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This final report was prepared by the Research Applications Division of 
Systems Research Laboratories, Inc., under NASA Contract NASI -16562, 
with Dr. G. Burton Northam as NASA Project Monitor. The research was 
conducted by Dr. Richard R. Antcliff. This report describes efforts 
performed during the period 12 March 1981 through 2 October 1986. 
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INTRODUCTION 


Physical probes have provided combustion researchers with temperature 
and pressure information for many years. However, interference with the 
flow properties and degradation in hostile environments limit the use of 
probes. Optical diagnostic techniques generally overcome these limita- 
tions. One of these techniques, coherent anti -Stokes Raman scattering 
(CARS), has the additional advantages of (1) high conversion efficiency, 
(2) a laser-like beam for high collection efficiency, fluorescence, and 
luminosity discrimination, and (3) high spatial and temporal resolution. 
These advantages must be weighed against the increased complexity of a 
CARS system and its limitations. These limitations include the presence 
of an interfering nonresonant background, the need for optical access, 

and the requirements for stability of the optical system. CARS has been 
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applied to many complex environments including plasmas, ’ internal com- 
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bustion engines, ’ shock tubes, * full-scale combustors, and a 

14 

supersonic combustor. 

This report describes a system which has been developed at NASA Langley 

Research Center for simultaneous measurement of temperature and nitrogen 

and oxygen number density. The system has been described, in part, pre- 
15 18 

viously as specific subsystems were developed. This system has 
been utilized to obtain a data base for comparison with computational - 
fluid-dynamics (CFD) calculations. These data have been included to 
illustrate the capabilities of the system. 

THEORETICAL BACKGROUND 

A CARS signal is generated as two laser beams at frequency oi-i and one 

laser beam at frequency 0)2 are mixed to produce a fourth laser-like beam 

at the CARS frequency = 2 uji - U 2 . The CARS signal is dependent upon 

the polarizability of the molecule. The polarizability can be written 
. 19 

as a power series 
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where E is the electric field strength and is the susceptibility. 
It is the third-order terms which describe CARS and other nonlinear 

Raman techniques. The intensity of the CARS signal, I3, has been 
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shown to be 
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where I-| is the intensity at frequency ; C is the speed of light; n3 

is the refractive index at 0)3; and L is the interaction length. This 

equation assumes that the required momentum phase matching has been 

satisfied and that the beams have the same focal diameter. The non- 
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linear susceptibility is a complex term which has been simplified for 
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normal CARS conditions * 


(3) _ N r . (3) 
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where N is the number density of the species probed; j is the popula- 
tion difference between the states studied; do/dn is the Raman cross 
section and assumed to be constant; is the frequency of a Raman 
active rotational or vibrational state; , is the Raman halfwidth; and 
xA ^ is a collection of susceptibility terms which are generally a 
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slowly varying function of frequency. The Xn ' term can become important 
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under resonance conditions, when near other rotational or vibrational 
transitions,^^ or when the species density is low.^^ 

From Eq. (3) it can be seen that as the difference between and 0)2 
approaches a Raman resonance, the signal increases dramatically. Nor- 
mally 0)^ is held fixed, and is scanned to generate a CARS spectrum. 

However, in turbulent environments this is impractical since the sample 
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changes rapidly. An alternative method makes use of a broadband dye 
laser for ^ 2 - Utilizing an array detector allows the complete spectrum 
to be acquired in a single 10-ns laser shot. 
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CARS SYSTEM 


The CARS system at NASA Langley employs a Nd:YAG laser, broadband dye 
lasers, nonplanar BOXCARS beam alignment, a reference signal generated 
in ambient air, and both photomultiplier- tube (PMT) and photodiode-array 
detection. Figure 1 is a schematic diagram of the optical layout and 
will be referenced in the following description. 

Laser System 

The main pump laser was a Quanta-Ray DCR-1 Nd:YA6 laser. This laser 
is normally operated at 10 Hz and emits 'v* 200 mJ in a 10-ns pulse at 
532 nm. The laser has been updated with new pump cavities and 0-ring 
seals. The 532-nm output was produced with a Quanta-Ray HG-2 
temperature-controlled frequency doubler. The temperature control is 
helpful in eliminating repeated tweaking of the crystal angle as the 
crystal heats up. The 532-nm output was separated from the residual 
1064-nm output by means of a Pell in-Broca prism. The residual output 
was directed through a second frequency doubler to obtain an additional 
532-nm laser beam. This portion of the 532-nm beam was used to pump the 
580-nm dye-laser oscillator. 

Two broadband (-v/ 200 cm"^ ) dye lasers were employed in this system. The 
oscillator of these dye lasers consists of a flowing dye cell oriented 
at Brewster's angle relative to the dye-laser optical axis. The oscil- 
lator was pumped slightly off axis in a cavity formed by a total and a 
30% reflector. A single amplifier stage was used for each dye laser. 
Several combinations of dye-laser concentrations and pump-laser energies 
were used in the oscillators and amplifiers to ascertain the optimum 
beam quality, power output, and stability. No advantage was found over 
using a simple 50/50 split of the laser energy between the oscillator 
and amplifier and the same dye concentration in both. 

The nitrogen-dye-laser (606.5 nm) concentration was -v 7.5 x 10'^ molar 
rhodamine 640 in methanol. The oxygen-dye-laser (508 nm) concentration 
was 5 X 10"^ molar rhodamine B in ethanol. The output frequency of 
these lasers can be adjusted by altering the concentration of the dye. 
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Figure 1. CARS Schematic. DC - doubling crystal, R - dichroic mirror, F - filter, 
A - polarization analyzer, T - telescope, l/2x - half-wave plate, 
l/4x - quarter-wave plate, CL - cylindrical lens. 




Signal Generation 
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A nonplanar BOXCARS ’ beam configuration was used to maximize the 
spatial resolution of the system while maintaining optimum signal- 
generation efficiency. Figure 2 shows the laser-beam configuration 
about the focus. The 532-nm "pump" laser was split into two equal beams 
and positioned vertically on the focusing lens. The two dye-laser beams 
were positioned vertically but placed to one side of the "pump" beams. 
The dye-laser beams should be slightly outside an imaginary circle on 
the lens having a diameter equal to the "pump"-beam separation and a 
midpoint at the center of the lens. With this configuration the signal 
beams were spatially separated from the input laser beams, and the input 
beams could be readily trapped. The crossing and focusing element was a 
plano-convex 50-cm lens. 

A small nitrogen jet was used to obtain the experimental interaction 
length. This jet, surounded by an argon sheath, was stepped through the 
interaction length and the signal recorded by measuring the signal to 
10% on either side of the peak. The effective signal -generation length 
of 1 mm was determined. The interaction volume was assumed to be a 
cylinder having a diameter of 250 ym. 

Referencing 

Measurements obtained from a stable environment (e.g., ambient air) 
showed considerable fluctuation in the signal intensity. This fluctua- 
tion was due to several factors including shot-to-shot pump-laser fluc- 
tuations (both in intensity and frequency), variations in the broadband 
dye laser output due to mode competition, and movement of the beams 
causing changes in the interaction length and in the detector response. 
This fluctuation presented a problem in making concentration measure- 
ments since such measurements are dependent on the signal intensity. To 
alleviate this problem a referening system was established to compensate 
for variation in the signal -generation efficiency. A 20% beam split- 
ter was placed immediately after the crossing and focusing lens. This 
splitter, placed at 45 deg., formed a second focal volume in ambient 
air. The nitrogen and oxygen signals generated in this reference volume 
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and the signals generated from the data (flame) volume were detected 
simultaneously. A flat glass plate was placed immediately after the 
splitter in the reference leg to compensate for a path difference in the 
sample leg caused by the splitter. The reference leg was optically 
adjusted to produce optimum correlation between the sample and reference 
signals. This was accomplished through use of a real-time program which 
displayed the correlation as the system was adjusted. Typical correla- 
tions were better than 10% for the nitrogen signals and better than 4% 
for the oxygen signals. 

Nitrogen Detection 

The complete nitrogen Q-branch spectrum was generated with a broadband 
dye laser. An intensified-sil icon-photodiode (ISPD) array was used in 
this system to capture a complete spectrum with each laser shot. The 
electronic scanning of this array of diodes was controlled by a mini- 
computer as described below. A 1-m monochromator was used in front of 
the detector to disperse the input beam and, thereby, reject stray laser 
light. 

A study was made to determine the experimental dynamic range of the 
1 8 

ISPD. Using calibrated sources the actual dynamic range of the 
instrument was found to be 100. Two major factors contributed to this 
small measured range. First, detector noise limited the lower measur- 
able limit to about 200 counts. Secondly, the spherical lens normally 
used to focus the CARS signal onto the detector was found to cause local 
saturation of the microchannel plate intensifier since the full detector 
height was not utilized. This dynamic range was insufficient for track- 
ing the large signal-intensity variations associated with turbulent com- 
bustion. 

Two modifications were made to the system to increase this dynamic 
range. First, the spherical lens at the entrance to the detector mono- 
chromator was replaced by a cylindrical lens. This lens produced a line 
image on the detector, thereby utilizing the full detector height. This 
improvement reduced the local saturation of the microchannel plate 
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intensifier and improved the linear dynamic range by approximately a 

factor of five. The second improvement was the implementation of an 
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external splitter arrangement as suggested by Goss. The nitrogen CARS 
signal was directed through a series of splitters to produce three sepa- 
rate signal beams having intensities of 1, 19, and 80% of the original 
signal. The monochromator entrance slit was rotated to the horizontal 
position and opened fully. The three signals were imaged side by side 
onto the slit plane. The signals followed the path of the monochromator 
and were imaged side by side on the detector. This allowed the three 
intensity signals to be recorded simultaneously by the diode array 
detector and increased the dynamic range by a factor of 100, which was 
sufficient for these studies. The reference nitrogen signal was also 
imaged onto the horizontal slit. Thus, four separate spectra were 
simultaneously recorded by the array detector. 

Oxygen Detection 

RCA 8575 Photomultipliers (PMTs) were used to detect the oxygen signals 
since only the integrated intensity of the signal was required. Mea- 
surements of the oxygen density are complicated by the presence of a 
nonresonant background signal at low oxygen densities. In this system 

the nonresonant background signal has been nearly eliminated by use of 
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polarization selection. The polarization characteristics of the reso- 
nant and nonresonant signals are different due to the electronic origin 
of the nonresonant signal. To take advantage of this difference, the 
polarization of the oxygen dye laser was rotated 60 deg. relative to the 
pump polarization. This was accomplished with a Karl Lambrecht 1/2-wave 
Fresnel Rhomb and a polarizing prism. A zero-order 1/4-wave retarder 
was also placed in the beam path to minimize elliptical polarization 
introduced by multiple optical elements. A polarizing prism was placed 
in the pump beam to improve the polarization characteristics of the 
beam. 

With the scheme discussed above, the polarization of the resonant signal 
is the same as that of the pump beams (Fig. 3). The nonresonant signal, 
however, was rotated 30 deg. from the pump beam. The nonresonant signal 
was blocked by placing an analyzer orthogonal to the nonresonant 
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required for rejection of nonresonant background. 


polarization axis. Some rejection of the resonant signal occurred in 
this configuration since the analyzer and resonant-polarization axis 
were separated by 60 deg. 

Bandpass filters were used to block stray light from the PMTs. Some 
initial complications were found due to the large bandpass of these fil- 
ters. A rotational hydrogen line (5 - 7) was passed by these broad 
rejection filters and detected by the oxygen PMTs. This signal was 
eliminated through the use of narrower bandpass rejection filters. 

System Calibration 

The accuracy and precision of the CARS measurements were demonstrated by 

use of a flat-flame burner. This burner had been studied previously and 
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found to be characterized by a well-behaved temperature profile. The 

burner, produced by McKenna Industries, Inc., was composed of a 7.94-cm- 

diam. sintered-stainless-steel plug housed in a stainless-steel jacket. 

The burner was cooled by a water coil within the plug. The burner was 

operated with the air flow held constant at 100 Ji/min and the hydrogen 

flow adjusted over the stoichiometric fraction range of 0.1 - 1.0. 

The water flow rate and temperature rise as well as the hydrogen and air 

flow rates were measured and recorded by the computer system. These 

experimental values were used to calculate the expected temperature and 
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species concentrations in the flow. 

Figure 4 shows the comparison of the CARS and calculated values for tem- 
perature, nitrogen number density, and oxygen number density. The CARS 
values shown represent the mean of 100 data shots. The oxygen measure- 
ments have a minimum detectable limit of 1 x 10^^ molec./cc due to the 
interaction of the nonresonant and resonant susceptibility which could 
not be completely eliminated because of incomplete polarization of the 
laser beams. 

The relationship between the intensity in the data legs and reference 
legs was also required for accurate calculation of species concentra- 
tion. This was accomplished by acquiring data with ambient air in both 
the data and the reference sample volumes. 
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COMPUTER CONTROL 


A complete description of the computer control of the Reticon detector 
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has been given previously; therefore, only a summary will be given 
here. The multichannel detector used in this system was an EG&G PAR 
1420 Reticon. This detector contains 1024 individual diode elements, 
the center 700 of which were preceded by a fiber-optic-coupled, 
proximity-focused, microchannel plate intensifier. An EG&G PAR 1218 
Detector Controller was used for electronic control of the scanning and 
data acquisition of the detector. A Reticon interface was specially 
built to communicate with the diode-array controller. This interface 
handles separate input and output data lines, signal drivers for long- 
distance operation, and the inverted logic required by the controller. 

Scanning of the Reticon was controlled by random-access memory (RAM) 
within the 1218 controller. This memory provided each individual detec- 
tor element (pixel) with predefined scanning instructions. These 
instructions were calculated by the MODCOMP and downloaded into the 
1218. This instruction set fast scans over the unintensified portion of 
the detector and provided normal analog-to-digital conversion over the 
center 700 pixels. 

An assembly-language program was written to initiate the Reticon scan- 
ning. At the beginning of each scan cycle, a signal was sent from the 
1218 to fire the laser. Three complete scans of the Reticon were accom- 
plished between laser shots. The first was the data-acquisition scan, 
and other two were used to clean the pixels of any remaining signal 
before the next laser firing. During these "clean scans" the data were 
transferred to magnetic tape for processing. 

Collection of photomultiplier data was also triggered by the laser pulse 
and, therefore, occurred simultaneously with the Reticon scans. A 
Kinetics 1500 CAMAC Crate was used to house the PMT collection electron- 
ics. The PMT signals were sent to a LeCroy 2249SG 12-Channel Analog-to- 
Ditigal Converter. This unit was triggered with a LeCroy 2323 Program- 
mable Dual Gate Generator. These signals were routed to the MODCOMP via 
an IEEE interface and stored on magnetic tape with the Reticon data. 
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An additional program was run simultaneously with the data-acquisition 
program to produce a real-time display of several system parameters. In 
particular, burner flow conditions and PMT output were displayed. The 
Reticon output was displayed on an oscilloscope. 

Data Reduction 


To obtain the temperature, the experimental nitrogen spectra were com- 
pared with a library of precalculated theoretical spectra. These spec- 
tra were calculated in 20-K steps and stored on disk. The experimental 
data were retrieved from tape and corrected with a background data scan 
and a detector response function. The computer then selected the larg- 
est unsaturated sample from the three nitrogen data signals. The data 
were compared with a library of precalculated theoretical spectra to 
obtain the experimental temperature. 

Once the temperature had been obtained, the concentrations of nitrogen 

and oxygen were extracted. The number density N of a species at tem- 
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perature t has been shown to be 


^t ^ '^300^^^t''b00^^^t,calc/^300,calc^^^^^ 

where N^qq is the number density at 300 K; is the experimental inte- 
grated intensity corrected for reference fluctuations; and I^qq is the 
integrated intensity at 300 K, obtained from a calibration run in room 
air. The calculated values were generated from the precalculated spec- 
tral library. 

The experimental nitrogen integrated intensity was obtained by summing 
the spectral channels selected by the fitting routine. Oxygen density 
calculations were simplified by the use of photomultiplier tubes which 
yield integrated intensity directly. A time of 2.5 s was required for 
analysis of each measurement on the MODCOMP computer. The arithmetic 
mean, standard deviation, and 95 % confidence level were calculated for 
each data set (100 laser shots) for comparison with CFD calculations. 
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CARS Signal Errors 


Errors must be considered which were associated with signal generation 
and signal detection. The most prevalent error associated with the CARS 
system is the fluctuation in the laser outputs. First, the YAG laser 
intensity varies considerably from shot to shot. Since this laser is 
used to pump the dye lasers, this variability is transferred to their 
output. Also, temperature changes in the doubling crystals effected a 
long-term change in the intensity output. A temperature-stabilization 
unit was installed on the primary doubling crystal to reduce this 
effect. In addition, referencing of the generation efficiency was used 
in this system to compensate for this type of fluctuation. The spectral 
character of the broadband dye lasers was dependent upon the fluores- 
cence and absorption of the dye used. However, because of the competi- 
tion among the many available modes, the output was very structured and 
variable. This type of fluctuation cannot be tracked adequately by 
referencing and constituted a major source of noise. 

In the current system, four laser beams were focused and crossed to form 
the sample volume. When sampling is conducted in the center of a turbu- 
lent environment, considerable beam steering may occur prior to the 
focus. This steering causes random misalignment of the beams and subse- 
quent random decrease in signal intensity. As the signal beams exit 
through the sample, beam steering may occur again. This steering may 
move the beam on the detection optics and alter the collection 
efficiency. 

The size and position of the sample volume was also critical to the 
interpretation of the data. In the burner studied, eddies of reactants 
and products form and move downstream. These eddies have a broad size 
range. If the sample volume is larger than the eddy size, an averaging 
of the temperature and densities will occur, and true point measurements 
will not be possible. In addition, if the sample volume is positioned 
such that it measures two different flows (i.e., fuel and coaxial air), 
a weighted averaging will occur. It should be noted that a tradeoff 
exists between sample-volume size and absolute signal intensity. 
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The nonresonant signal can be a source of error. The nitrogen analysis 
software assumed that the nonresonant signal was relatively low and, 
therefore, handled in the data reduction. The oxygen system, because of 
the relatively low oxygen density, utilized polarization rejection of 
the nonresonant background. This technique requires excellent polariza- 
tion characteristics in the lasers and optics for achievement of com- 
plete rejection. In this system distortions in the polarization char- 
acteristics of the laser beams were introduced by the large number of 
optical elements through which the laser beams had to pass. Therefore, 
the possible incomplete rejection of nonresonant background limited the 
minimum oxygen level detected. 

The narrow filter used before the oxygen PMTs to reject an interfering 
hydrogen line (see above) may have distorted the lineshape of the oxygen 
signal. This distortion could have caused an error in the oxygen- 
density calculation since it is dependent upon a known lineshape. The 
errors associated with the Reticon detector used in the nitrogen leg 
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have been well documented. These errors were found to be due mainly 
to saturation of the image intensified on the front end of the photo- 
diode array. 


CARS DATA-REDUCTION ERRORS 

The theoretical description of the CARS process was continuously being 
refined to fit the experimental data more accurately. These improve- 
ments are generally not needed for such low-resolution studies. How- 
ever, they may introduce finite errors and, therefore, are discussed 
here. In this system the nonresonant-background term was assumed to be 
a simple function of constituent densities which were adiabatically cal- 
culated from the experimental temperature. The nonresonant values of 
Rado^^ were used. These values have been disputed by Hall.^^ Recently, 

more complex fits which include the nonresonant background as a variable 
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have also been attempted with success. This system also assumed that 

the Raman linewidth follows a simple 1/T function. Rahn has shown that 
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this assumption is an oversimplification. The theoretical spectra 

must be convolved with the system laser and detector linewidth func- 
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tions. The convolutions of Yuratich have been used for this system. 
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Recently, Teets and Katolica introduced an improved convolution 
which more accurately takes into account the nonresonant interference. 
The library interval (20 deg.) also imposes a lower limit on the reduced 
data. 


SUBSONIC DIFFUSION BURNER 

The burner used in the present studies consisted of two coaxial tubes. 
The inner tube, which carried the fuel, had a 4.7-mm inner diameter and 
a 6.4-mm outer diameter. This tube was tapered at the exit plane to a 
lip thickness of 0.51 mm. The outer tube, which carrried air, had an 
inner diameter of 24.5 mm. The length of the burner was 45 cm. Typi- 
cal flow rates were 100 m/s in the air tube and 15 m/s in the fuel tube. 
The fuel consisted of a mixture of 80% hydrogen and 20% nitrogen. The 
addition of nitrogen to the fuel was necessitated by the use of nitrogen 
as the temperature probe species. 

Figure 5 displays a spark shadowgraph of this burner under the above- 
stated flow conditions. This obviously very turbulent flow provided a 
good test of the dynamic range of the CARS system. Fast photography of 
this flow reveals the formation of large eddys which break up further 
downstream. 

The subsonic burner was mounted on a motorized X-Z translational system. 
The Y-axis (along the laser beams) was adjustable but not motorized. 

This system effectively allows the sample volume to move through various 
positions of the flame without alignment difficulties. 

The normal acquisition sequence began with a calibration run in ambient 
air. After the burner had been started and stabilized, a background run 
was made. Each data point was the accumulation of 100 laser shots. 
During the background run the dye-laser beams were blocked to allow only 
stray light to enter the detectors. The burner was then positioned by 
means of the motorized translational stages, and a CARS data point was 
obtained. One additional background point was obtained for approxi- 
mately every five data points. 
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Figure 5. Shadowgraph of combustion produced by coaxial diffusion burner. 






An axial profile and two radial profiles were obtained above the coaxial 
burner. The axial profile began as close to the jet exit as physically 
possible and extended 35 fuel-jet diameters downstream (Fig. 6). The 
radial surveys were located at 5.5 and 22.1 fuel -jet diameters down- 
stream (Figs. 7 and 8). These radial profiles began at the center line 
and extended 5 fuel -jet diameters to the outer portion of the flame. 

COMPUTATIONAL FLUID DYNAMICS 

The major application of the CARS results to date has been in the vali- 
dation of computational -fluid -dynamics (CFD) computer-model ing codes. 
Considerable effort has been placed on generating codes for predicting 
complex combustion environments. CARS measurements provide experimental 
data which can be used to refine the numerous adjustable parameters 

found in these codes. Initially, two separate codes were used with the 
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data. The first, a modification of a CHARNAL code, solves parabolized 
Navier-Stokes equations using a marching finite-difference algorithm and 
includes a two-equation turbulence model (k - c). Combustion was mod- 
eled by hydrogen-oxygen equilibrium chemistry, with the nitrogen being 
treated as inert. Input to the code consisted of velocity, fuel mass 
fraction, turbulent kinetic energy, static temperature, and static pres- 
sure at the initial station located 0.01 fuel-jet diameters downstream 
of the exit plane. Initial velocity and turbulence intensity profiles 
were obtained with a hot-wire anemometer using cold-air flows. Initial 
profiles of static temperature and pressure were assumed to be uniform 
over both jets at 300 K and 1 atm. Equilibrium chemistry assumes that 
the local mixture reacts instantaneously to an equilibrium mixture of 
products; this assumption is typically used with low flow rates. 

The second computational method used was the TEACH code, which is a 

modified version of the REACT2D^^ code. The code solves the elliptical 

time-averaged Navier-Stokes equations using the k-e turbulence model for 

two-dimensional reacting flows. The eddy-breakup combustion model of 
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Magnussen and Hjertager was used; this model relates the combustion of 
fuel to the breakup and mixing of turbulent eddies and the mean concen- 
tration of reacting species. A solid wall located at 4 jet diameters was 
used to approximate the air boundary. The law of the wall velocity 
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INE PROFILE 




Figure 6. CARS data and CFD calculation along centerline. 


RADIAL PROFILE 
5.5 Dia. downstream 



Figure 7. CARS data and CFD calculations from radial profile 5.5 jet diameters 
downstream. 


RADIAL PROFILE 
22.1 Dia. downstream 
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Figure 8. CARS measurements and CFD calculations from radial profile 22.1 jet diameters 
downstream. 



profile was used. At the outlet, axial gradients of all variables 
except velocity were assumed to vanish. The TEACH code iterated until 
the nondimensional residual was minimized. 

In general, the correlation between the CARS data and the data generated 
by the CFD codes was very good. The major differences were in the flame- 
spreading rate and the peak temperature. These deviations are readily 
explainable by the assumptions built into the CFD codes. Two factors 
may contribute to this deviation— an overprediction of the heat release 
and an overprediction of the mixing rate. The heat release may be 
over-predicted for two reasons. First, the equilibrium model may not 
adequately model the real system which has finite reaction rates. Sec- 
ondly, these codes do not include a complete coupling model between 
turbulence and the chemical reaction. The mixing rate may have been 
overpredicted for three additional reasons. First, numerical diffusion 
can cause numerical smearing; this effect was suspected to be minor 
since the flow was aligned along the mesh and refinements to the grid 
had only small effects upon the results. Secondly, the turbulence model 
itself, even for nonreacting flows, overpredicted the mixing of the jet. 
Thirdly, as the jet reacts, a low-density "buffer region" forms in the 
shear layer; the turbulence model did not account for the dynamic 
effects of this variable density. 

The CHARNAL code appears to furnish more accurate predictions of the 
mixing as indicated by the relatively close correlation between the 
experimental and calculated data in the axial profile and the lower 
predicted spreading rate in the radial profiles. However, the peak tem- 
peratures calculated by the TEACH code more closely match the experimen- 
tal data. The lower peak temperature may be the result of the use of 
the eddy-breakup model. This model will weakly link the chemistry and 
turbulence, thereby preventing the system from obtaining the "equilib- 
rium" temperature value. 


SUMMARY 

A nonintrusive diagnostic system has been developed around a coherent 
anti -Stokes Raman spectroscopy instrument. This system is capable of 
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measuring temperature and nitrogen and oxygen number densities in turbu- 
lent combustion environments. High temporal and spatial resolution, 
luminosity discrimination, and a large dynamic range are important 
advantages of this system. Turbulent subsonic flow has been investi- 
gated with this system. 
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